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ABSTRACT 

We use high-resolution A-body simulations to follow the formation and evolution of tidal 
streams associated to dwarf spheroidal galaxies (dSphs). The dSph models are embedded 
in dark matter (DM) haloes with either a centrally-divergent ‘cusp’, or an homogeneous- 
density ‘core’. In agreement with previous studies, we find that as tides strip the galaxy the 
evolution of the half-light radius and the averaged velocity dispersion follows well-defined 
tracks that are mainly controlled by the amount of mass lost. Crucially, the evolutionary tracks 
behave differently depending on the shape of the DM profile: at a fixed remnant mass, dSphs 
embedded in cored haloes have larger sizes and higher velocity dispersions than their cuspy 
counterparts. The divergent evolution is particularly pronounced in galaxies whose stellar 
component is strongly segregated within their DM halo and becomes more disparate as the 
remnant mass decreases. Our analysis indicates that the DM profile plays an important role in 
defining the internal dynamics of tidal streams. We find that stellar streams associated to cored 
DM models have velocity dispersions that lie systematically above their cuspy counterparts. 
Our results suggest that the dynamics of streams with known dSph progenitors may provide 
strong constraints on the distribution of DM on the smallest galactic scales. 

Key words: galaxies: dwarf, kinematics and dynamics, evolution, Local Group - cosmology: 
dark matter - methods: numerical 


1 INTRODUCTION 

It is widely believed that the microscopic properties of dark mat¬ 
ter (DM) particles are reflected upon the distribution of DM on the 
smallest galactic scales (e.g. Tremaine & Gunn 1979; Gilmore et al. 
2007). Being the smallest and most dark matter-dominated galax¬ 
ies in the known Universe, dwarf spheroidal galaxies (dSphs) place 
strong constraints on the DM particle mass (e.g. Boyarsky et al. 
2009; Maccio et al. 2012; Lovell et al. 2014) as well as on its cross 
section (Burkert 2000; Loeb & Weiner 2011; Vogelsberger et al. 
2012, 2014; Rocha et al. 2013). Furthermore, with luminosities in 
the range of 10 4 — 1O 8 L 0 (McConnachie 2012), these galaxies sam¬ 
ple the faintest end of the galaxy luminosity function, providing a 
testbed for theoretical models that aim to understand the dynamical 
interplay between baryons and the DM haloes (e.g. Navarro et al. 
1996; Pontzen & Governato 2012, 2014; Penarrubia et al. 2012; 
Brooks & Zolotov 2014; Di Cintio et al. 2014; Nipoti & Binney 
2015). 

However, measuring the distribution of DM from the inter¬ 
nal kinematics of dSphs has proven to be challenging. Their low 
surface brightness plus their relatively large heliocentric distances 
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limit the number of stellar tracers for which kinematic informa¬ 
tion can be obtained and also reduce the phase-space information 
to quantities projected along the line of sight. In addition, as fore¬ 
ground contaminants tend to dominate in number in the outskirts 
of these galaxies most spectroscopic surveys suffer from strong 
sampling biases that severely complicate the dynamical analysis 
of these objects. Given the number and graveness of these short¬ 
comings it is perhaps not surprising that different methods that aim 
to measure the mass profile of Milky Way dSphs yield contradic¬ 
tory results. For example, whereas some studies find evidence for 
constant-density DM cores (Battaglia et al. 2008; Walker & Penar¬ 
rubia 2011; Amorisco & Evans 2012), other methods often applied 
to the same data suggest centrally-divergent ‘cusps’ (Richardson & 
Fairbairn 2014) or remain inconclusive (Breddels & Helmi 2013; 
Strigari et al. 2014). 

Recent numerical studies have shown that the dynamical prop¬ 
erties of tidal streams are sensitive to the internal mass distribu¬ 
tion of the progenitor system (Kiipper et al. 2008, 2012; Penarrubia 
et al. 2006, 2010; Eyre & Binney 2011). Using high-resolution N- 
body simulations we show in this paper that the distribution of dark 
matter in dSphs plays an important role in defining the orbits of 
tidally stripped stars, and that the dynamics of stellar streams pro¬ 
vide constraints on the dark matter halo profile of dSphs which are 
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Figure 1. Evolutionary tracks for cuspy and cored dSph galaxies. All quan¬ 
tities are measured at apocentre (see text). The dashed lines are fits of Equa¬ 
tion 2. The half-light radius r^, the luminosity average velocity dispersion 
Oh within r/* and the total stellar mass m* are normalised to the initial value 
and are plotted as a function of the current mass mh/m^o enclosed within 
r h . The tracks are independent of the orbit, slightly dependent on the ini¬ 
tial stellar segregation r*/a and significantly dependent on the shape of the 
mass profile (cuspy/cored) of the dSph galaxy. 


Table 1 . Empirical Fit Parameters to the Tidal Evolutionary Tracks 
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complementary to those inferred from the internal kinematics of 
the progenitor galaxy. The Letter is organised as follows: Section 
2 describes the Af-body dSph and analytical host galaxy models as 
well as the numerical setup. In Section 3 we analyse how cuspy and 
cored dSphs evolve under tidal stripping and introduce the evolu¬ 
tionary tracks. In Section 3 we present our results regarding the dif¬ 
ferences of tidal stream originating from cuspy and cored models. 
Section 4 describes how these differences can be used to constrain 
the DM profile of dSphs with associated tidal streams. 


2 NUMERICAL METHOD 

The host galaxy is modelled as a Milky Way-like static potential 
composed of a spherical NFW halo (Navarro et al. 1997) with 
virial radius r ylv = 258 kpc, mass M ylv = 1O 12 M 0 and concentra¬ 
tion c = 12 (Klypin et al. 2002); a spherical Bulge (Hernquist 


1990) with scale radius a^ = 1.2 kpc and mass = 1.3 • 1O 1O M 0 , 
and an axis-symmetric disk (Miyamoto & Nagai 1975) with ra¬ 
dial and vertical scale lengths a j = 3.5kpc, bj = 0.3kpc and mass 
Md = 7.5 • 10 10 M o . 

To generate A-body models of dSphs in dynamical equilib¬ 
rium we draw 2 x 10 7 particles from a Dehnen (1993) distribution 
function. These models follow an analytical density profile 


M (3 —y) a 
47tr7(r + a) 4- ^ 


( 1 ) 


where y = — lim r ^o dlogp /dlogr defines the slope of the inner den¬ 
sity profile. In this work we consider two values, namely y = 1 
(cusp) and y = 0 (core), and set a = 1 kpc. Our dSphs models are 
built under the assumption that stars behave as massless tracers of 
the underlying DM potential. To generate a stellar component em¬ 
bedded in the DM halo in equilibrium we follow the approach de¬ 
scribed by Bullock & Johnston (2005), where each individual DM 
particle is assigned a probability, or mass-to-light ratio, of repre¬ 
senting a ’star’ in such a way that the overall stellar density pro¬ 
file follows a Plummer sphere with segregation r*/a~ 7^/(1.3a), 
where r^ and r * are the half-light and core radii of the Plum¬ 
mer model, respectively. To facilitate a comparison between cored 
and cuspy models the dSph parameters are chosen such that for 
a given spatial segregation both cored and cuspy DM haloes have 
the same mass enclosed within the half-light radius nth = m(< r^) 
at t = 0. Since the mass profile of a Dehnen model is m(< r ) = 
M[r/(r + a)] 3 ~Y, imposing m^ cusp = m/^ core requires changing the 
total mass of the dwarf models by M cusp /M core = r^j(r^ + a). Fix¬ 
ing M cusp = 10 8 M o yields M core = 8.66 • 1O 8 M 0 for r*/a = 0.1, 
and M CO re = 4.83 • 1O 8 M 0 for r*/a = 0.2. Our models are purely 
gravitational, thus they are scale free, and physical units are just 
given for clarity. 

We use the particle-mesh code Superbox (see Fellhauer et al. 
2000), which samples the density in three grids with different reso¬ 
lution and performs a leapfrog integration to solve the equation of 
motion. In our models each grid consists of 128 3 cubic cells. Grids 
1 and 2 move with the centre of density of the progenitor, while 
grid 3 is centred on the host galaxy. Grid 1 resolves the core of the 
progenitor with a resolution of 2a/126 « 16pc, while grids 2 and 
3 have resolutions of 20a/ 126 « 160pc and 400a/126 « 3.2kpc, 
respectively. We use a fixed time step of 1 Myr and store snapshots 
every 200 Myrs. All models are run in isolation for 14 Gyrs in order 
to test the stability of our numerical set-up. 

The dSph models are placed on highly eccentric orbits with 
a fixed apocentre distance r a = 150kpc and pericentres r p = 
5,7.5,10,15 and 20kpc. These orbits are confined to the plane per¬ 
pendicular to the symmetry axis of the host disc in order to facil¬ 
itate the analysis of the associated tidal streams. As tides strip the 
dSph models we follow the evolution of the half-light radius r^, the 
luminosity-averaged velocity dispersion o\ and the stellar mass m * 
as a function of the mass m/* within r^. All quantities reported in this 
paper are non-projected. Note that can be in principle measured 
with relatively high accuracy in isolated systems (see e.g. Walker 
et al. 2009; Wolf et al. 2010), but the calculation of <7 requires a 
luminosity-averaged integral over the entire dSph, which in case 
of the dSph undergoing tidal stripping may be difficult to estimate 
owing to the presence of transient ‘extra-tidal’ features. To compute 
the half-light radius and the stellar mass m*, we fit a Plummer 
profile to the DM particles weighted by the corresponding mass-to- 
light ratio. To avoid contamination with ‘extra-tidal’ particles we 
follow an iterative approach, where only 3/4 of the particles within 
a radius of 5 kpc to the dSph are fitted at a time. Then using the fit- 
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Figure 2. Snapshots of the evolution of tidal streams associated to cuspy (top row) and cored (bottom row) models. Both models have r a = 150kpc,r p = 
20kpc,a = lkpc and r*/a = 0.2. One orbital period equals roughly 2.4Gyrs. The masses have been rescaled to ensure that the luminosity-averaged velocity 
dispersion within r\ % is the same for both progenitors at t = 0. Small boxes in each panel show a detail view of the progenitor. The rightmost column shows the 
projection of tidal stream and progenitor on the (R,z) plane, where R is the cylindrical distance from the galactic centre and z is the distance above the orbital 
plane of the progenitor. DM particles are colour-coded according to their probability of tagging stars. See ArXiv ancillary files for animations 1 . 


ted Plummer parameters we derive the radius that contains 3/4 of 
the stellar mass and repeat the fit again. These steps are re-iterated 
until convergence is reached. The stellar velocity dispersion and 
mass nth are computed from the velocities of the particles located 
within the half-light radius. 

Because the internal dynamical time of the dSph models are 
much shorter than the period of the orbit, the progenitor reaches 
an equilibrium configuration shortly after each pericentre passage. 
To decrease stochastic fluctuations parameters are averaged over 7 
snapshots chosen symmetrically around each apocentre. 


3 TIDAL EVOLUTIONARY TRACKS 

Previous studies (e.g. Penarrubia et al. 2008,2010) have shown that 
tidal effects tend to lower the velocity dispersion and increase the 
half-light radii of stars embedded within Navarro et al. (1997) ha¬ 
los, and that the structural changes in the stellar component of these 
galaxies follow well-defined evolutionary tracks that are mainly 
controlled by how much (and not how) the mass m^ enclosed within 
the half-light radius varies. Fig. 1 shows that this result also holds 
for galaxies embedded in cored DM haloes. However, while the 
tidal tracks associated to cuspy DM haloes scarcely depend on the 
spatial segregation of the stellar component within the DM halo, as 
shown by Penarrubia et al. (2008), those associated to cored DM 
haloes appear to be sensitive to the initial value of r±/a. In particu¬ 
lar, we find that the half-light radius increases to a larger extent the 
more deeply segregated the stars are within the DM core. For the 
remainder of this paper, it is useful to parametrize the evolutionary 
tracks using the empirical formula of Penarrubia et al. (2008) 

2^jc^ 

/(*)= (l +x ) a WherC X = m h/ m h,0 ■ (2) 

The best-fitting parameters are listed in Table 1 and the resulting 
tracks are plotted with dotted and dashed lines in Fig. 1. Fig. 1 also 
shows other points of interest. First, note that at a fixed remnant 
mass nth/mhfi cored progenitors have larger half-light radii and 


colder velocity dispersions than their cuspy counterparts. This dif¬ 
ference can be easily understood in terms of the dynamical time of 
the progenitor galaxies, which scales as ^yn ^ (Up) -1 / 2 . Because 
Pcusp Pcore at r a we find that stars embedded in DM cusps 
react adiabatically to external perturbations. In contrast, stars em¬ 
bedded in cored DM haloes tend to absorb a relatively large amount 
of energy during tidal encounters, which leads to an overall expan¬ 
sion and a low kinematic temperature once dynamical equilibrium 
is reached. 

A second interesting point can be gleaned from the lower 
panel of Fig. 1, which shows the fraction of stars lost to tides ver¬ 
sus the remnant mass. We find that at a given remnant mass cored 
dSph models lost less stars than their cuspy counterparts. This ap¬ 
parently counter-intuitive result arises from the dichotomy between 
the shapes of the stellar and DM density profiles. Whereas stars in 
in cored DM halo models tend to be associated to the most-bound 
DM particles, in cuspy DM haloes a large fraction of particles with 
high binding energy are located within the density cusp and have 
a low probability of tagging a star, i.e. a high mass-to-light ratio. 
E.g. for cuspy models with r^ja — 0.2, one finds that 90% of the 
stellar particles are associated to the 14.5% most-bound DM parti¬ 
cles, while in cored models this fraction falls to a mere 5%. Note, 
however, that because dSphs embedded in cored DM haloes are in 
general less resilient to tidal stripping, these models tend to evolve 
at a faster pace along the evolutionary tracks than those embedded 
in cuspy haloes (e.g. Penarrubia et al. 2010). 


4 STREAM EVOLUTION 

In this section we compare the properties of stellar streams associ¬ 
ated to cored and cuspy progenitors. As expected, Fig. 2 shows that 
tides are strongest at pericenter and trigger recurring episodes of 
mass loss at every pericentric passage. Comparison between the 

1 dSph_stream_brightness.avi: colour-coded by surface brightness 
dSph_stream_dispersion.avi: colour-coded by velocity dispersion 
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Figure 3. Width w and velocity dispersion a of the streams normalized to 
the half-light radius and luminosity-averaged velocity dispersion of the pro¬ 
genitor dSph, respectively. We show data at a fixed snapshot (4 th apocentre) 
of the simulation shown in Fig. 2. We distinguish between dSph models 
embedded in cuspy (open symbols) and cored (filled symbols) DM haloes. 
Symbols are colour-coded according to the pericentric passage at which 
particles became tidally unbound. We mark the cylindrical distance R to the 
galactic centre of segments within the range 20kpc ^ R ^ lOOkpc. Notice 
the clear kinematic offset of the streams depending on the DM halo profile. 


upper and lower panel reveals that tidal mass loss is more pro¬ 
nounced in dSphs embedded in cored DM haloes. Indeed, although 
not shown here the cored dSph is fully disrupted during the 5th 
encounter, whereas the cuspy dSph still retains a bound remnant. 

At each snapshot we quantify the mass The width 

w 2 = (z 2 ) — (z) 2 and the velocity dispersion g z of particles moving 
along the tidal tails are measured in the direction per pendicu lar to 
the orbital plane at fixed cylindrical distances R = \/ x 2 + y 2 from 
the host galaxy centre. These quantities are then normalized by the 
half-light radius and the velocity dispersion of the stream progeni¬ 
tor as given by Equation 2 for the corresponding value of m^/m^ q. 
To reduce noise we apply the same averaging is it was done with 
the tracks. Fig. 3 shows that tidal streams become very narrow and 
kinematically hot at small distances from the host and get wider 
and colder with increasing R , as also shown in Fig. 2, as one would 
expect from Liouville’s theorem. 

In order to understand the internal evolution of the stream 
properties it is useful to distinguish between particles lost at dif¬ 
ferent pericentric passages. We define a wrap as the ensemble of 
particles stripped at the same pericentre encounter. Each wrap is 
then divided into segments according to their cylindrical distance R 
to the galactic centre. Wraps that formed during recent pericentric 
passages tend to be found at close distances from the progenitor 
galaxy, as it takes a few orbital periods for these particles to pop¬ 
ulate the entire orbit. In contrast stars stripped at earlier times can 
in some cases wrap around the host galaxy entirely. Comparison 
between different wraps in Fig. 3 shows that streams associated 
to dSphs embedded in cored DM haloes have similar widths but 
higher velocity dispersions than their cuspy counterparts. This is so 
regardless of the galactocentric radii at which these quantities are 
measured. Also, wraps that form during early encounters tend to 
have higher velocity dispersions than those formed recently. This 
trend arises because the velocity dispersion of the dwarf galaxy de¬ 
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Figure 4. Time evolution of the width and velocity dispersion ratios be¬ 
tween streams associated to dSphs embedded in cored and cuspy dark mat¬ 
ter haloes of the simulation shown in Fig. 2. Symbols denote particles lost at 
different pericentric passages (wraps), whereas colours distinguish between 
different snapshots. The width and velocity dispersion ratios roughly follow 
the ratio between the tidal evolutionary tracks of cuspy and cored dSphs as 
given by Equation 2, which we plot with a dotted line for ease of reference. 


creases at a faster pace (see Fig. 1) than the velocity dispersion 
of the associated stream. Hence, tidal streams appear kinematically 
hotter with respect to to the progenitor’s velocity dispersion as they 
become older 2 . Fig. 3 shows that both the width and the velocity 
dispersion of tidal streams are highly dependent on the progenitor’s 
orbit, the number of orbital periods completed, and the galactocen¬ 
tric radius where these quantities are measured. 

Interestingly, Fig. 4 shows that plotting the ratios <7 C ore/(Tcusp 
and Wcore/^cusp at fixed segments collapses the degeneracy be¬ 
tween these effects into a single relation, which roughly follows 
the path defined by the tidal evolutionary tracks of the progeni¬ 
tor galaxies, Equation 2, evaluated for the corresponding values of 
m h/ m h ,0 f° r the cuspy and cored models (shown with dotted lines 
for ease of reference). This result suggests that the evolutionary 
path of the stream in this plane is mainly driven by the disparate 
structural evolution of galaxies embedded in cored and cuspy DM 
haloes as a result of tidal mass loss (see Fig. 1). As shown in 
Fig. 4, there is an intrinsic offset between the various properties of 
the stream which can be understood using simple dynamical argu¬ 
ments. If we assume, as in Amorisco (2014), that the stream width 
is proportional to the tidal radius of the progenitor at pericentre 
r t oc M 1 / 3 , and that the velocity dispersion of the stream scales as 
G z \/m(< r t )/r t , one finds that for the orbits plotted in Fig. 4 
the width and velocity dispersion ratios are r^ CO re/^,cusp — 1.7 and 
t^core/o^cusp ~ 1.7, respectively (see section 2). Fig. 4 shows that 
these simple analytical estimates, marked with a black cross for 
ease of reference, faithfully describe the numerical results. It is 
straightforward to show that for models with smaller stellar segre¬ 
gation r*/a with the same values of r h and G^ this offset increases. 
This is in agreement with our models for r*/<z = 0.1 (not shown 
here). 

2 Stream ‘age’ is typically defined as the time elapsed since the particles 
become tidally unbound. 
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5 CONCLUSIONS 

We have used A-body simulations to follow the evolution of dSphs 
as driven by galactic tides. Our models adopt cosmologically- 
motivated mass profiles where dSphs are embedded in extended 
dark matter haloes with density profiles that have either a centrally - 
divergent ‘cusp’ or an homogeneous-density ‘core’. As in previ¬ 
ous studies, we find that the half-light radius r^, velocity dispersion 
Oj t and stellar mass m * of dSphs acted on by tides follow well- 
defined evolutionary tracks that are mainly controlled by the rem¬ 
nant mass nth/nth ,o enclosed within the half-light radius. Crucially, 
these tracks also appear to be remarkably sensitive to the adopted 
shape of the halo mass profile (cusp/core), and also to a minor ex¬ 
tent to the degree of segregation of stars within the dark matter 
halo. 

Our models indicate that the differences in the evolutionary 
paths are reflected upon the internal dynamics of the associated 
tidal streams. In particular, after normalizing the stream width and 
velocity dispersion by the half-light radius and internal dispersion 
of the remnant progenitor, respectively, we find that tidal streams 
of dSphs embedded in cored dark matter models tend to be system¬ 
atically hotter than their cuspy counterparts. The kinematic offset 
arises during pericentric encounters and becomes more pronounced 
as the age of the stream increases. 

The results presented in this Letter call for a simple method 
for inferring the distribution of dark matter in dSphs with associ¬ 
ated stellar streams. First, one would derive by measuring 

the stellar mass in the stellar stream and in the remnant galaxy. For 
a given dark matter halo profile and a given initial stellar segrega¬ 
tion this observation will fix the position of the galaxy on the tidal 
tracks shown in Fig 1, and thus the properties of the progenitor 
system prior to tidal stripping. Models that reproduce the observed 
rh and Oh lead to M core > M cusp . Having set the initial conditions 
of the dwarf galaxy, the next step corresponds to modelling the 
orbit and the mass evolution of the progenitor system simultane¬ 
ously in a given host potential. This can be done with the aid of 
self-consistent N-body models (e.g. Bonaca et al. 2014) or, alterna¬ 
tively, with efficient semi-analytical algorithms based on forward 
modelling of stellar streams (e.g. Varghese et al. 2011; Sanders & 
Binney 2013; Amorisco 2014; Gibbons et al. 2014; Fardal et al. 
2014; Price-Whelan et al. 2014). Our results suggest that the ve¬ 
locity dispersion of stellar streams associated to cuspy halo models 
will lie systematically below that of the cored models, which in turn 
can be used to constrain the mass profile of the progenitor galaxy 
by comparing the velocity dispersion of the stream against that of 
the remnant galaxy. 

To date the only dwarf spheroidal with a known tidal stream 
corresponds to the Sagittarius dSph (Ibata et al. 1997, see Be¬ 
lokurov et al. 2014 for a recent data compilation.) As this dSph 
is currently located close to pericentre, we can’t apply directly our 
tidal evolutionary tracks, as these have been computed for dSphs 
that have reached equilibrium after each pericentre passage. The 
use of TV-body simulations will be required to study perturbations 
from the equilibrium configuration. An application of the proposed 
method to the existing kinematic data set will be presented in a 
separate contribution. 
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